Computing the Ball Size of Frequency Permutations 

under Chebyshev Distance 



Min-Zheng Shieh 
Department of Computer Science 
National Chiao Tung University 
1001 University Road 
Hsinchu City, Taiwan 
Email: mzhsieh@csie.nctu.edu.tw 



Shi-Chun Tsai 
Department of Computer Science 
National Chiao Tung University 
1001 University Road 
Hsinchu City, Taiwan 
Email: sctsai@csie.nctu.edu.tw 



o 

(N 



> 

o 



X 



Abstract — Let S,^ be the set of all permutations over the 

multiset {1,...,!,... , m, . T. , rri} where n — mX. A frequency 
permutation array (FPA) of minimum distance d is a subset 
of 5*;^ in which every two elements have distance at least d. 
FPAs have many applications related to error correcting codes. 
In coding theory, the Gilbert- Varshamov bound and the sphere- 
packing bound are derived from the size of balls of certain 
radii. We propose two efficient algorithms that compute the ball 
size of frequency permutations under Chebyshev distance. Both 
methods extend previous known results. The first one runs in 

O (^Cax)^ ''^'^'^ iogn) time and O ((^J^^)^) space. The second one 

runs in O {{^ax)i x^)j) time and O {{^^^)) space. For small 
constants A and d, both are efficient in time and use constant 
storage space. 



I. Introduction 

Frequency permutation arrays (FPAs) of frequency A 
and length n are sets of permutations over the multiset 



{1, . . . , 1, . . . , TO, . . . , m} where n = mX. In particular, FPAs 
of frequency A = 1 are called permutation arrays (PAs). 
FPA is a special case of Slepian's codes [10.1 for permutation 
modulation. Many applications of FPAs in various areas, such 
as power line communication (see [|8], ifTsl . lHJ^ and fTsI), 
multi-level flash memories (see lO and fH, ifTTl ) and computer 
security (see lfT2l ). are found recently. In many applications, 
we use FPAs as error correcting codes. It is well-known that 
the capability against errors of a code is mainly determined 
by its minimum distance. Similar to traditional codes, the 
minimum distance of an FPA C is d under metric 5 if 
minp.7rec;p7^7r ^{p, tt) — d. A (A, n, d)-FPA under some metric 
(5 is a FPA of frequency A and length n which has minimum 
distance at least d under S. A (A, n, d)-FPA is often considered 
to be better if it has larger cardinality. When we evaluate the 
quality of certain design of (A, n, d)-FPA, we often compare 
it with the maximum cardinality of {X,n, d)-FPAs. Generally 
speaking, computing the maximum size of codes of certain 
parameters is hard. In coding theory, the Gilbert- Varshamov 
bound and the sphere-packing bound are famous lower and 
upper bounds on the code size, respectively. They are derived 
from the size of balls of certain radii. We focus on the 
efficiency of computing the ball size in this paper 



Shieh and Tsai lfT2l showed that the cardinality of d-radius 
balls can be obtained by computing the permanent of a matrix. 
It is #P-complete to compute the permanent of a general 
matrix. However, the matrix for estimating the ball size of 
FPAs has a special structure. Kl0ve |5| used the property and 
proposed a method to efficiently compute the cardinality of 
balls of radius 1 by solving recurrence relations. Meanwhile, 
based on Schwartz's result f9\, one can compute the size of 
d-radius balls in O (^{^^) '^^^logrt^ time for A = 1. 

In this paper, we give two algorithms to compute the cardi- 
nality of a d-radius ball under Chebyshev distance for d > 1 



and A > 1. The first one runs in O 



O 



f2dX\ 



(( 



dX . 



2.376 , 



log n ) time and 



<1dX\ ldX-\-X\ n 
\dx)\ A ) X 



J I space, and the second one in O 

time and O (^{^^^x)^ space. These algorithms are generaliza- 
tion of Schwartz's result [9j. They are efficient in time and 
space when d and A are small fixed constants. 

The rest of the paper is organized as follows. In section 
2, we define some notations. In section 3, we introduce a 
recursive algorithm EnumV for enumerating permutations in 
a d-radius ball. Then, based on the property of EnumV, we 
give two methods to obtain the ball size efficiently in section 
4. We compare previous results with ours in section 5. Then, 
we conclude this paper briefly. 

II. Notations 



We set n — m\ throughout the work unless stated other- 
wise. For positive integers a and b with a <b, [a] represents 
the set {1, . . . , a} and [a, 6] represents {a, a + 1, . . . , 6 — 1, &}. 
For convenience, (— oo, &] represents the set of integers at 
most h. The Chebyshev distance of two fc-dimensional vectors 
X and y is defined as d,„ax(a;,y) = maxigf^j \xi - yi\, 
where Xi and are the i-th entries of x and y respectively. 
For permutations x and y, they are said to be d-close to 
each other under metric 6{-,-) if 6{x,y) < d. We use 
to denote the set of all frequency permutations with each 
symbol appearing A times. The identity frequency permutation 



e in 5„ 

e is Ci 



is (1, . . . , 1, . . . , TO, . . . , rri), i.e., the i-th entry of 
= [ j] . A partial frequency permutation can be 



derived from a frequency permutation in S^^ with some entries 
replaced with *. The symbol * does not contribute to the 
distance. I.e., the distance between two fc-dimensional partial 
frequency permutations, x and y, is defined as dynaxix, y) — 
maxjg[i.] J. .^^ y .^^ \xi — yi\. 

Under Chebyshev distance, a ball of radius r centered at tt 
is defined as B{r,Tr) — {p : doo{p,'^) < r}. We can obtain 
B{r, tt) from B{r, e) by rearranging the indices of the entries, 
therefore \B{r, tt)\ ^ \B{r, e)| for tt e S^. Let Voc>iX, n, d) be 
the size of a ball of radius d in under Chebyshev distance. 

III. Enumerate rf-CLOSE permutations 



In this section, we give a recursive algorithm to enumerate 
all frequency permutations in B{d,e). First, we investigate e 
closely. 
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Observe that symbol fc appears at the (fcA — A + l)-th, . . . , 
(fcA)-th positions in e. Therefore, tt is d-close to e if and only 
if TTfeA-A+i, . • . , TTfcA G [k — d, k + d]. (Note that we simply 
ignore the non-positive values when fc < d.) In other words, 
rfmax(e, tt) < d if and only if symbol fc only appear in the 
(fcA — dA — A + l)-th, . . . , (fcA + dA)-th positions of vr. This 
observation leads us to define the shift operator ©. 

Definition 1. For an integer set S and an integer z, define 

S ® z ^ {s + z : s e S}. 

Fact 1. Suppose that tt is d-close to e, then -Ki — fc implies 

i e [-d\ + 1, dA + A] ® (fcA - A). 

The above fact is useful to capture the frequency per- 
mutations that are d-close to e. Note that the set S = 
[— dA + 1, dA + A] is independent of fc. We give a recursive 
algorithm ENUMVA,„,d in figure[T| It enumerates all frequency 
permutations tt G S^^ in B{d,e). It is a depth-first-search 
style algorithm. Basically, it first tries to put I's into A proper 
vacant positions of tt. Then, it tries to put 2's, . . . , m's 
into the partial frequency permutations recursively. According 
to fact [U symbol fc is assigned to positions of indices in 
[— dA + 1, dA + A] © (fcA — A), and these positions are said to 
be valid for fc. 

ENUMVA,n,d takes an integer fc and a subset P of [— dA + 
1, dA + A] as its input, and ENUMVA.ri.d uses an {n + 2dA)- 
dimensional vector tt as a global variable. For convenience, we 
extend the index set of tt to [— dA+1, n+dA] and every entry of 
TT is initialized to 0, which indicates that the entry is vacant. We 
use P to trace the indices of valid vacant positions for symbol 
fc, and the set of such positions is exactly P © (fcA — A). 

We call ENUMVA,n,ci(l, [dA+A]) to enumerate all frequency 
permutations which are d-close to e. During the enumeration, 
EnumVa, n, ti(fc, -P) assigns symbol fc into some A positions, 
indexed hy X (B (fcA — A), of a partial frequency permutation, 
then it recursively invokes ENUMVA,,i,ci(fc + 1,X' ® (-A)) U 
[dA + 1, dA + A]), where X and X' form a partition of P and 
\X\ = A. After the recursive call is done, EnvmY \^n.d{k, P) 



ENUMVA,„,rf(fc,P) 

1. if fc < TO then 



2. for each partition {X, X') of P with \X\ = A do 

3. if X' n (-oo, -dA + A] = then 
// Make sure it is a proper partition 

4. forieX® (fcA - A) do 

5. TTi <— fc; // Assign fc to the i-th position 

6. y ^ (X' © (-A)) U [dA + l,dA + A]; 

7. ENUMVA,„,d(A: + 

8. for i e X © (fcA - A) do 

9. TTi <~ 0; // Reset tt^ to be vacant 



10. else 

11. if P = [dA + A] then output (tti, . . . , 7r„); 

Fig. 1: ENUMVA,„,d(fc,P) 

reset positions indexed by X © (fcA — A) as vacant. Then, 
it repeats to search another choice of A positions until all 
possible combinations of A positions are investigated. For 
fc = TO + 1, ENUMVA,„,d(fc, P) outputs TT if P = [dA + A]. 
Given n = Ato, fc is initialized to 1 and P to [dA + A], we 
have the following claims. 

Claim 1. In each of the recursive call of EnumV a. ri,d, in line 
6 we have max(y) = dA + A. 

Proof: By induction, it is clear for fc = 1. Suppose 
max(P) = dA + A. Since Y ^ {X' (S (-A)) U [dA + 1, dA + A] 
and max(X') < ^-^ + we have max(y) = dA + A. ■ 

Claim 2. In line 6, for each fc G [to + 1] and each i G 
y © ((fc + 1)A - A), we have tt., = 0. 

Proof: We prove it by induction on fc. It is clear for k — I. 
Assume the claim is true up to k < m + I, i.e., for each 
i G P©(fcA— A), TTi = 0. Now, consider the following scenario, 
ENVM'Vx,n.d{k,P) invokes EnumV x,n.d{k + l,Y). 

Since Y = {X' ® (-A)) U [dA + 1, dA + A], we have Y © 
((fc + l)A-A) = (X'©(fcA-A))U[fcA + dA + l,fcA + dA + A]. 
While X' C P and [fcA + dA + 1, fcA + dA + A] are new vacant 
positions, it is clear TTi = in these entries. ■ 

Claim 3. In each recursive call of EmMW \ n,d{k, P), P must 
be a subset of [—dA+1, dA+A] of cardinality dA+A. It implies, 
- k\ < dfor i e P ® (fcA - A). 

Proof: We prove it by induction on fc. For fc = 1, P is 
[dA + A], and the claim is obvious. Assume the claim is true 
up to fc, and ENUMVA,„,d(fc, -P) invokes ENUMVA,„,d(fc + 
l,y). Thus y = {X' ©'(-A)) U [dA + l,dA + A]. Due to 
the constraint on X' in line 3 and the induction hypothesis, 
we have X' © (-A) C [-dA + 1, dA]. We conclude that Y C 
[-dA+l,dA]U [dA + 1, dA + A] and |y| = \X'\ + X = dA + A. 
Since [-dA + 1, dA + A] © (fcA) = [fcA - dA + 1, fcA + dA + A], 
we know e has values from [fc — d + 1, fc + d + 1] in these 
positions. I.e., [e^ — (fc+ 1)| < d for i G y © (fcA). Hence, the 
claim is true. ■ 

Claim 4. At the beginning of the invocation of 
ENUMVA,„.d(fc, P), « G P © (fcA - A) implies i > 0. 



Proof: It is clear for k = 1. Observe that min(F) > 
min(P)-A. Since, (min(P) - A)©(fcA) = min(P)©(fcA-A), 
the claim holds for fc > 1. ■ 

Claim 5. For k £ [m], when E^VMW \^n.d{k,P) invoke 
ENUMVA.n,!i(fc + 1;^) ^^e 7, there are exactly A entries 
of TT equal i for i G [fc — 1]. 

Proof: It is implied by lines 4 and 5. ■ 

Lemma 1. At the beginning of the execution of 
ENUMVA.„.rf(fc, P), P ® (fcA - A) = {i : i > A TT, = A i G 
[-fiA + l,dA + A]® (fcA-A)}. 

Proof: The lemma holds by claims |2] |3] and |4] ■ 

Let Trbe the output of ENUMVA,n,d(l, [c?A+A]). These facts 
ensure that tt G 5*,^ and (imax(e, tt) < d. 

Lemma 2. For fc G [m + 1], let tj. be a partial frequency 
permutation d-close to e and with each symbol 1, . . . , fc — 1 
appearing exactly A times in Tk- For i = let 
T^i = {Tk)i, if {Tk)i 7^ *■ Let TTi = for the other i. Every fre- 
quency permutation p in S^^, generated by EnumVa, n,ci(fc, P), 
satisfies the following conditions: 

1) p is consistent with Tk over the entries with symbols 
l,...,fc-L 

2) p is d-close to e. 

Proof: We prove it by reverse induction. First, we con- 
sider the case fc m + 1. By claim |5] every symbol appears 
exactly A times in r,„+i and tt. By claim |2] and claim ID 
we know that all of tti , . . . , 7r„ are nonzero if and only if 
P = [dX + A]. There are two possible cases: 

» P — [d\ + A]: By claim [51 (tti, . . . , 7r„) is the only 
frequency permutation in satisfying both conditions. 

• P ^ [d\ + A]: Then there is some i G [n] such that 
TTi = 0. But there is no frequency permutation p in 
with Pi =111 + 1 > m. This means vr is not well assigned. 

Note that ENUMVA,,i.(i(fc, -P) outputs tt only if P = [(iA + A], 
otherwise there is no output. Hence, the claim is true for 
k = m + 1. Assume the claim is true down to fc + 1. For fc, by 
lemma [T] P©(fcA — A) is exactly the set of all positions which 
are vacant and valid for fc. In order to enumerate permutations 
satisfying the second condition, we must assign fc's into the 
valid positions. Therefore, we just need to try all possible 
choices of A-element subset X©(fcA— A) C P©(fcA — A). Line 
3 of EnvmW \,n,d{k, P) ensures that X is properly selected. 
Then symbol fc is assigned to X © (fcA — A) and we have a new 
partial frequency permutation Tk+i- By induction hypothesis, 
it is clear that the generated frequency permutations satisfy 
both conditions. ■ 
We have the following theorem as an immediate result of 
lemma |2] 

Theorem 1. ENUMVA,n.(i(l, [dA© A]) enumerates exactly all 
d-close frequency permutations to e in S^^. 

Proof: Let p be a d-close frequency permutation to e in 
S^. There is a corresponding sequence of partial frequency 



permutations ti, . . . ,Tm+i d-close to e. Note that ti has * 
in all of its entries and Tm+i = p- By lemma |2l p will be 
enumerated eventually. ■ 

IV. Computing Ko (A, ri,d) 

The number of elements generated by ENUMVA,n,(i(l, [dA© 
A]) is clearly Voo{X, n,d). However, the enumeration is not 
efficient, since FooCA, is usually a very large number In 
this section, we give two efficient implementations to com- 
pute V^tx)(A, n, d). Especially, 14o(A, n, d) can be computed in 
polynomial time for constant d and A. 

From the algorithm EnumV, we see that whether 
ENUMVA,„.d(fc,-P) invokes ENUMVA,„.d(fc + 1,^) or not 
depends only on fc, P and Y . During the execution of 
ENUMVA,„,rf(l, [d\ + A]), ENUMVA,n,d(fc, -P) is invoked re- 
cursively only when [d\ + l,d\ + \] C P C [-dX + l,dX + X], 
due to line 6. Therefore, we can construct a directed acyclic 
graph Gx,n,d = {Vg,Eg) where 

. Vg = {(fc,f/) : fc G [m + !],[/ C [-dX + l,dX] and 
\U\ =dX}. 

. {{k,U),{k + 1,V)) G Eg if and only if 

ENUMVA,„,d(fc, U [dX + l,dX © A]) invokes 
ENUMVA,„,d(fc + l,VU[dX + 1, dX © A]). 




Fig. 2: Graph 6*2,6,1 



For example, figure |2] shows the structure of G2.6,i- 

Voo{X, n,d) equals the number of invocations of 
ENUMVA,„.<i("i © l,[dX © A]). With this observation, 
Voo{X, n, d) also equals the number of paths from (1, [dX]) to 
(m © 1, [dX]) in G\,n,d- By the definition of G\,n,d, it is a 
directed acyclic graph. The number of paths from one vertex 
to another in a directed acyclic graph can be computed in 
0{\V\ + \E\), where \V\ = (to©!)^^'^^) and \E\ = 0{\V\^). 
So Voo(A, n,d) can be calculated in polynomial time with 
respect to n if A and d are constants. 

The computation actually can be done in O(logn) for 
constant A and d. Define Hxd = {Vh,Eh) where Vh = {P : 
|P| = dX and P C [-dX © 1, dX]} and (P, P') G Eh if and 
only if there is some fc G [m] such that ENUMVA,n.d(fc, -P U 
[dX © l,dA © A]) invokes ENUMVA,„,d(fc © 1,P' U [dX © 
l,dX © A]). Note that \Vh\ = {^^^j- kgure [3] shows H2.1 
as an example. 




the number of edges has an 



Fig. 3: Graph H2,i 



Theorem 2. Voo{X,n, d) can be computed in O(logn) time 
for constant d and A. 

Proof: Observe that the value of A; e [to] is independent 
of the invocation of ENUMVA,„,d(fc + 1, P' U [dA + 1, dA + A]) 
by ENUMVA,„,d(fc, FU[dA + l,dA + A]) , where P and P' C 
[—dX + IjdA] with |P| = \P'\ = dX. Therefore, the number 
of paths of length to from [dX] to itself in H\ d is equal to 
the number of paths from (1, [dX]) to (to + 1, [dX]) in G\,n,d- 
Let Vh = {wi, ■ • ■ , v\yjj\}, where vi — [1, dA]. The number 
of paths of length to from vi to vi is the first entry of 
the first column of the TO-th power of Ah, where Ah is 
the adjacency matrix of Hx d- Since TO-th power can be 
computed in O i^f ^(^^_!^)^ logTO^, where 0{f{x)) is the time 
cost of multiplymg two x y. x matrices. It is well-know that 
f{x) — 0{x'^^^'^^) by the Coppersmith- Winograd algorithm. 
With constants A and d, V\^n,d can be found in 0(log n) time. 

■ 

However, the space to store the adjacency matrix Ah, 
^ij^dx) )' ^® large to execute the 0(logn)-time 
algorithm. For example, by setting d — Z and A = 3, we need 
at least ( g^) w 2.36 x 10^ entries to store This makes the 
constant factor extraordinarily large in the proof of theorem |2] 
Hence, we provide an alternative implementation which runs in 
O {(^^^) ■ C^Y^) • time and 0((^/a)) space. This allows 
us to compute more efficiently for the cases with smaller to 
and larger d and A. For example, (^^^) ■ C^^^^) • to w 10^ 
for d = 3, A = 3, and TO = 100. To achieve the 0{{^^^))- 
space complexity, we do not store the adjacency matrix Ah 
in the memory. Since Ah is the adjacency matrix of H\^d, 
fory = {yi,...,y\VH\) and y' = AhV = {y[, . . . ,y'^y^^), we 
have y'^ — J2{ij)GEH Vi' H^i^ce, if enumerating all edges in 
Eh takes S space and T time, then we can compute Any in 
0(|Vfr| + S) space and 0{T) time for any iVf/j -dimension y. 



Lemma 3. \Eh\ < \Vh\- C^^x^) and Eh can be enumerated 
in 0{dX) space and 0{\Eh\) time. 



has 



Proof: For P GVh such that |Pn(-oo, -dA-A] \ = r, P 
out-going edges, since every partition {X, X') 



fd\+\-r\ 
A-r 



of P U [dX + 1, dA + A] satisfies the condition in line 3 if and 
only if Pn(— oo, — dA — A] C X, i.e., every choice of (A — r)- 
element subset of oo, — dA — A] will invoke a recursive 



call. Since C^Y-rl ^ 
upper bound IV^hIC'Y^)- 

To enumerate all A-element subsets of a (dA + A) -element 
set, we need 0(dA + A) = O(dA) space and 0((''^+^)) time. 
Since we can recycle the space, the enumeration of edges in 
Eh can be done in O(dA) space and 0{\Eh\) time. ■ 
Now, we give the alternative implementation. 

Theorem 3. Vqo{X, n,d) can be computed in 

o {i'^x) ■ i'T) ■ ™) ^p"^^- 

Proof: Let x = (1, 0, . . . , 0)^. Since A'gx is the first 
column of A^, Voo{X, n,d) is the first entry of A^^x. The 
alternative evaluates . . . , A^a; iteratively. Instead of 

storing the whole adjacency matrix, it only uses two \Vh\- 
dimension vectors y and y' for storing A^a; and the interme- 
diate result of A'^f^^x, respectively. Initially, y = x and i = 0. 
We compute Any by the algorithm described in lemma |3] and 
using y' to store the intermediate result. Then, we copy the 
result of Any back to y. After to repetitions, we have y = 
A'^x, and the first entry of y is Voo(X, n,d). Therefore, the 
space complexity can bereducedtoO(|Vi/|+dA) = 0{{^^^^)). 

The running time is to • 0{\Eh\) = O (^{^^^) ■ ('^^+^) • to^ 

V. Comparison with previous results 

In this section, we compare our results with previous ones. 
Shieh and Tsai 1121 showed that (A!)'"Voc.(A, n, d) equals 
the permanent of 0-1 matrix y^C^ ".'*) 



^^{\ n,d).^ where 



a. 



{X.nA) 
i,3 



1 if and only if | [j] — [^] | < d. For example. 



^(2.8,1) 



A naive approach to evaluate the permanent of an n-by-n 
matrix takes 0{n\) time. In practice, 9((A!)'"ycxj(A, n, d)) 
time is still required when using a backtracking algorithm. 
It is clear that both of our method are much more faster 

Kl0ve fSl solved the recurrence of Voo (A, n, d), and he gave 
the value of T4o(A, Ato, d) for A G [10], to G [20], and d = 
1. Schwartz ||9l gave an algorithm which can be applied for 
computing V'oo(l, n, d). In this paper, we provide solutions to 
computing VooiX^n^d) for A > 1 and d > 1, which is not 
contained in their works. We list the values of 14o(A, Ato, d) 
for A > 1, TO e [20], d > 1, and dA < 10. 

Voo(2, 2, 2) = 1 
Voo(2, 4, 2) = 6 
Voo(2, 6, 2) = 90 
Voz,(2, 8, 2) ^ 786 
Voo(2, 10, 2) = 6139 
Vcxj(2, 12, 2) ^ 54073 
Voo (2, 14, 2) — 477228 
Voo(2, 16, 2) = 4113864 
Voo (2, 18, 2) = 35579076 
Voo (2, 20, 2) — 308945881 
Voo(2, 22, 2) — 2679325561 
Voo (2, 24,2) —23222971098 
Voc(2, 26, 2) = 201351085146 
Voc(2,28,2) = 1745886520422 



: 15137227297027 
: 131243141767393 
: 1137923361184848 
. 9866167034815440 
: 85542686564024352 
: 741681846818742097 



: 9747523 

: 154700569 

: 2502207156 

: 40043708244 

: 632349938520 

: 9986116318524 

. 158192179607364 

: 2509767675626581 

: 39796612230719845 

: 630688880128338378 

: 9994168619297530758 

. 158396161513685960664 

. 2510580301930785916566 

: 39792149406721332018414 



: 90 

^ 2520 

^ 113400 

= 3540600 

= 88610850 

= 2044242426 

= 47806940971 

= 1196081134201 

= 30647443460124 

= 784921116539484 

= 19899840884886720 

= 500019936693729120 

= 12551808236761063440 

= 315694279415609776404 

= 7955400980632212027852 

— 200622722060793477132937 
= 5057787000067792980984649 

— 127452627155747602225756890 



: 90 

^ 2520 

= 113400 

— 7484400 

= 361859400 

= 14091630840 

= 489147860970 

= 16420511188146 

= 563209318269379 

= 20416518083009593 

= 758713036253909844 

= 28351365170599079604 

= 1054143198114097909680 

= 38864351069181445164480 

= 1423417411123883479886400 

= 52064892889568503574209920 

= 1906534315066176639758670480 

= 69931615009402042606373019S04 

^ 1 

. 20 

. 1680 

= 61340 

= 1886431 

— 69496201 

= 2568223000 

= 91712960320 

= 3290467596440 

= 118724053748417 

= 4276273204804217 

= 153904262366842444 

= 5541519231941145440 

= 199545071017172522244 

= 7184755645113714298863 

= 258691998154725997048673 

= 9314545233907934721851472 

= 335381528796576643131475840 

= 12075785123501322139824319056 

= 434802491356562053648077727185 

: 1 
: 20 
: 1680 

= 369600 

— 41480880 

— 3422150780 

= 276888204387 

= 25512718688405 

= 2418264595619240 

= 225661997838758560 

= 20649533952628896000 

= 1889648253594082624960 

= 173699198403114756474600 

= 16001577154624484682748453 

= 1472965856766989578006355117 

= 135481185586476496195656612044 

= 12459839493182349378716705969200 

= 1146141579672729885487800599057600 

= 105440511941055519854115528116882480 

= 9699923367172090411762252385134967844 

: 1 
: 70 

= 34650 

= 5562130 

= 708212251 

= 114774147001 

= 18679465660540 

= 2906167849870600 

= 454904037056013460 

= 71729455730285511001 

= 11285129375761977675001 

= 1773699532985462649188410 

= 278931562239767189408085850 

= 43869015908453746845566145990 

= 6898693708786029238293860809251 

= 1084865341390442288732669957148001 

= 170605963060816377946936433265175680 

^ 26829411396875692269491197638918648400 

^ 4219165662049303123773116859323196816720 

= 663502408038018748448058464247159216890001 



: 298227062281 
(5, 30, 2) = 218838390759073 
(5,35,2) ^ 161446400503248672 
(5,40,2) — 112632613848657302400 

. 79169699996993643966432 
■ 56151546386557366024202177 
^ 39717291593245217794329362081 
. 28058660061656964336359435570604 
: 19835819533825566529982592591911412 
^ 14024417724324420598672399947721245804 
(5,75,2) = 9914206081036463014882722168252570938889 
(5,80,2) = 7008596284293402975749309111124669929079521 

. 4954676885097638926007640423100194180529855296 
: 3502659589845301193905028251874353899223998638208 
: 2476160267409321946445662150301548547825713614803904 
(5, 100, 2) = 1750492069977099993617695861204414333904857504132837761 



Conclusion 

We extend Schwartz's result Q to solve the ball size 
of frequency permutations under Chebyshev distance. We 
also offer another efficient algorithm for the cases of larger 
frequency, larger minimum distance and smaller symbol set. 
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^ 252 

■ 756756 

i 549676764 



